Granular convection and the Brazil nut effect in reduced gravity 
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We present laboratory experiments of a vertically vibrated granular medium consisting of 1 mm 
diameter glass beads with embedded 8 mm diameter intruder glass beads. The experiments were 
performed in the laboratory as well as in a parabolic flight under reduced-gravity conditions (on 
Martian and Lunar gravity levels). We measured the mean rise velocity of the large glass beads 
and present its dependence on the fill height of the sample containers, the excitation acceleration, 
and the ambient gravity level. We find that the rise velocity scales in the same manner for all three 
gravity regimes and roughly linearly with gravity. 

PACS numbers: 81.05.Rm 44.25.+f 45.70. Mg 81.70.Ha 



The vertical segregation of particle sizes, also known 
as the Brazil nut problem [T] , has long been known but is 
very complex in its nature. A main reason for the com- 
plexity of the problem is that there exist different driving 
mechanisms for different experimental parameters, like, 
e.g., container shape and size, excitation acceleration and 
frequency, among others, and that these mechanisms are 
not mutually independent but can overlap and take place 
at the same time. One widespread driving mechanism 
is granular convection as shown by Knight et al. [21 [3] , 
who studied the convective motion of glass beads in long 
cylindrical Pyrex and Lucite containers by using dyed 
tracer particles as well as magnetic resonance imaging. 
With this mechanism, a collective motion of the medium 
can transport larger particles to the top, which come to 
rest there if the downward flow zone is too small to be 
entered by these particles. This effect has been studied 
experimentally [2 HTT| and theoretically ^121115] , while the 
focus of the article at hand is on the extrapolation of the 
Brazil nut problem to reduced gravity conditions found 
on small Solar System bodies. The Brazil nut effect has 
been, for example, made responsible for observed sur- 
face structures on small asteroids HH [T^ . Since many 
small bodies in the Solar System possess a granular sur- 
face [18] there are also technical aspects to this problem, 
like, e.g., the handling of materials on the Moon, mate- 
rial sampling on upcoming asteroid landing missions, or 
on larger scales even future asteroid mining (e.g., consid- 
ering density segregation). For all of these problems, it is 
desirable to understand the scaling of the granular flow 
and, thus, the segregation timescale with the ambient ac- 
celeration Jamb, which can be as small as 10~^(7Eaith on 
a small asteroid. 

We studied the Brazil nut effect of glass-bead samples 
in a transparent polycarbonate container of 110 mm in- 
ner diameter and 100 mm height. This container was 
firmly mounted on an infeed slide of a linear rail to verti- 
cally shake it in the direction of its cylindrical axis (Fig. 
fl] left). The shaking profile in all experiments had an 
oscillation amplitude of ±10 mm and was defined by an 
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FIG. 1: Left: Sketch of the experimental setup. The sample 
container is attached to an infeed slide, which can be moved 
vertically on a linear rail. Right: Acceleration profile of the 
infeed slide and the attached experiment container for an ex- 
periment with ^amb = SMars- 



approximated square-wave function of the acceleration 
(see Fig. [l] right). Depending on the individual experi- 
mental run, a pre-defined acceleration level was set to the 
control software of the linear stage and the acceleration 
was internally regulated. We additionally attached an ex- 
ternal acceleration sensor to the infeed slide to measure 
the acceleration level and found a small overshoot before 
reaching a quasi-constant acceleration plateau at the pre- 
defined value (Fig. IT right). We will use the maximum 
acceleration amplitude of the periodic accelerometer sig- 
nal as indicated by the gray-shaded areas for the exci- 
tation acceleration Ocxc- One should keep in mind that 
the vibration frequency was not kept constant with this 
setup, but changed from 7.0 to 11.1 Hz (for Earth- level 
gravitational acceleration), 5.1 to 5.7 Hz (Mars) and 3.4 
to 3.8 Hz (Moon) for the the different experimental runs. 



We used soda-lime glass beads of 1 mm diameter for 
the bulk material and 8 mm diameter for the intruders. 
The test chambers were initially filled to a height of a 
few millimeters with the small particles, then seven in- 
truder beads were fixed in this granular bed in a hexag- 
onal shape with one bead in the center, before the con- 
tainer was filled up with 1 mm beads to a nominal height 
of 60 mm in most cases. We do not expect an influence 
on the rise time due to multiple intruders, as confirmed 
by Mobius et al. [8]. The experiments were performed 
under normal air pressure and we therefore cannot ex- 
clude gas effects. However, due to the chosen particle 
size, mass density of the bulk material, and an equal den- 
sity between bulk and intruder particles, we expect that 
the particle rise time is only slightly affected [5]. The 
upper surface of the sample was observed with a digital 
camera with a precise internal timer to measure the rise 
time of the glass bead intruders, which we took with an 
accuracy of 1 s. A horizontal movement of the surface 
layer and some heaping was visible in the camera images. 
Moreover we observed that large beads, which appear at 
the surface, reemerge into the medium near the side walls 
and reappear after a time, which is always less than twice 
the rise time, so that we conclude to be in a convective 
regime. 

Some experiments were performed on ground and, 
thus, under Earth gravity (gEarth = 9.81 ms~^) and we 
measured the rise time (rise velocity) as a function of 
the fill height and the excitation acceleration, as will be 
presented below. The main focus of this work was, how- 
ever, on low-gravity experiments, which were performed 
on board the A300 Zero-G aircraft flying a maneuver to 
achieve Martian (.gMars — 3.71 ms^^) and Lunar grav- 
ity (gMoon ~ 1-62 ms^^) levels. In these experiments, 
the excitation of the granular sample was started when 
the desired gravity level was attained and the excitation 
was stopped shortly before the transition to the following 
hyper-gravity maneuver. This yielded a time between 20 
and 30 s per parabola and if the intruders did not ap- 
pear in the first parabola, the experiment was continued 
in the subsequent parabola. The granular medium was 
perfectly in rest when the excitation was stopped so that 
the experimental runtime can simply be cumulated. In 
few cases, the excitation was only stopped within the 
hyper-gravity phase, but at the abrupt increase of the 
gravity level, the granular medium always came to a com- 
plete rest, which was clearly visible on the camera images. 
In these cases, we chose the time at which the granular 
medium stopped moving at the end of the experimental 
runtime of that particular parabola. In each of the four 
fiights, in which 12 Mars- and 12 Moon-level parabolas 
were flown, we used several experiment containers, which 
were exchanged after all seven intruders had appeared. 

In a first set of experiments on Earth, we kept the am- 
plitude of the excitation acceleration constant at a nor- 
malized value of r = 1.72. Here, F = aoxc/ffamb is the ra- 
tio of the excitation acceleration and the ambient gravita- 
tional acceleration (in this case with gamb = ffEarth)- We 
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FIG. 2: The mean rise velocity as a function of the fill 
height of the experiment containers (red circles, left ordinate) . 
Also plotted is a cumulative distribution of fill heights in the 
parabolic flight experiments (green squares and blue triangles, 
right ordinate), which shows that the range of fill heights (58 
to 66 mm) was close to the maximum rise velocity. 



varied the fill height from 20 to 80 mm and placed the in- 
truder beads always at the bottom of the container. The 
mean rise velocity (i.e., the fill height minus the intruder 
diameter, divided by the rise time) is shown in Fig. [2] as 
a function of fill height. One experiment with seven in- 
truders is represented by one red circle with the standard 
deviation shown by the error bars. We found a linear in- 
crease of the rise velocity for fill heights up to 60 mm. 
This is consistent with the observation of Garcimartin 
et al. [B], who found the flow velocity (downward flow 
near the wall) in their experiments to scale linearly with 
the number of layers at aspect ratios (container height 
divided by container diameter) from 0.18 to 0.45 but at 
higher excitation frequencies of 110 Hz. For larger fill 
heights, we found a decrease of the rise velocity, which 
has been studied earlier by [Knight et al.| and [Grossman 
[31 [13] . These authors found an exponential decay of the 
rise velocity with increasing embedding depth for much 
larger aspect ratios than ours. For aspect ratios less than 
unity, they found that the rise velocity decreases with in- 
creasing fill height (but with an unknown relation), so 
that our results are also consistent with these findings. 
We are not aware of any publication showing an abrupt 
change in the rise velocity as presented in Fig. [5] The 
fill height for the Earth-gravity experiments was precise 
to within one millimeter, but due to different filling pro- 
cedures in the low-gravity parabolic fiight experiments, 
the latter show a small variation. Thus, we also plotted 
in Fig. [2] a cumulative distribution of fill heights of the 
parabolic-flight experiments, which are between 58 and 
66 mm and, thus, in a regime with maximal rise velocity. 
Figure [3] shows the rise velocity of the terrestrial and 
parabolic-flight experiments, in which the fill height was 
either 60 mm (Earth) or 58 to 66 mm (Mars and Moon) 
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FIG. 3: Rise velocity for experiments at terrestrial (red cir- 
cles), Martian (green squares), and Lunar gravity values (blue 
triangles) as a function of the excitation acceleration, here 
represented by F — 1 . The error bars show the standard devi- 
ation in one experimental run with seven intruders. The solid 
and dashed lines represent power laws with an exponent of 
1.3 as described in the text. 



as a function of the excitation acceleration F — 1. For 
the terrestrial experiments, where gamb = 5Earth, we find 
that the data are well reproduced by a power law of the 
form 



(X (F - 1) 



1.30 



(1) 



thus the choice of F — 1 as the horizontal axis. We did 
not use the data point with the highest F — 1 value for 
the fit as it clearly deviates from the other data for an 
unknown reason. A vanishing velocity at F < 1 is con- 
sistent with most earlier works [SHS]. To compare our 
results to the data of Knight et al. [3] , we translated our 
rise times into tap numbers, thus multiplied it with the 
oscillation frequency. Considering their data for 1 mm 
glass beads in Figs. 3a and 4a and their Eq. 1, we can 
see a rough quantitative agreement for an initial depth 
of z = 60 mm, but we can get a perfect match (slope 
and absolute value) for z = 30 mm. A possible reason 
to choose an initial depth smaller than our fill height is 
their smaller container diameter, which is smaller than 
ours by a factor of two. We could not find an agree- 
ment with other publications, which were however either 
sparsely described pXJj or used a 2D setup [II |9] . Avail- 
able data for the convection velocity at the walls, which 
show an exponent in Eq. [I] between 1 and 2 [SJ [TT] , can- 
not directly be compared to our data as Hejmady et al. 
[9] showed that the exponent for the near-wall velocity 
is close to unity, while the rise velocity (as derived from 
their rise times) possesses a much steeper dependence on 
excitation acceleration. 

The excitation acceleration and the rise velocity in the 
parabolic-fiight experiments was measured the same way 
as for the terrestrial experiments. However, since the 



FIG. 4: Rise velocity for F = 1.9 as a function of the ambient 
gravity level (?amb for the three gravity regimes studied. 



ambient gravitational acceleration (i.e.. Lunar or Mar- 
tian gravity) is not really constant (because the aircraft 
encounters turbulence) , we have to apply a correction to 
the rise-velocity data. Even though the mean accelera- 
tion over a full parabola comes very close to the nominal 
Lunar or Martian gravity level, an oscillation of typically 
iO.OS^Earth around this value has a non-negligible impact 
on the rise velocity. If the aircraft acceleration increases 
(decreases) over the nominal value, the relative excita- 
tion F gets smaller (higher). The corrected rise velocities 
|19) are shown as green squares (Martian experiments) 
and blue triangles (Lunar experiments) in Fig. pj again 
representing mean values for one container with seven in- 
truders. The scatter of these data is slightly larger than 
for the terrestrial experiments - most likely still due to 
variations in the ambient gravity level - but a clear trend 
of increasing rise velocity with increasing excitation level 
is evident. A power-law fit to the two low-gravity data 
sets yields exponents of 1.47 (Mars) and 1.17 (Moon), re- 
spectively, which are close enough to the terrestrial value 
of 1.30, which is based on a much wider range of excita- 
tion accelerations. The green and blue solid lines in Fig. 
[3] represent power laws with a slope of 1.30. The dashed 
lines represent the expected velocities for Martian and 
Lunar gravity levels if the scaling with gravity is linear. 
The Martian data are consistent with this linear extrap- 
olation, while the Lunar data significantly deviate from 
the linear trend. 

To derive the relation between the rise velocity and am- 
bient gravity level, we chose a reference point of F = 1.9 
for the normalized excitation acceleration (roughly the 
median acceleration of all experiments). For this nor- 
malization, each data point was shifted according to Eq. 
[TJand the highest F — 1 value for the Earth experiments, 
which was also not used for the fit, was again neglected. 
The resulting mean rise velocity of all experiments in 
one gravity regime is shown as one data point in Fig. HI 
(same symbols as in Fig. an for the three ambient gravi- 



tational accelerations. The error bars represent standard 
deviations of the mean values. One can see that the rise 
time is roughly in a linear relation to the ambient gravity 
level, as indicated by the solid line. The deviation of the 
rise time for the Lunar experiments is, however, statisti- 
cally significant so that we cannot make a firm statement 
on the functional behavior. A perfect fit can actually be 
achieved with an exponential function, shown by the dot- 
ted line, with the drawback that it gives an implausible 
non-zero rise velocity for a vanishing ambient acceler- 
ation, if gravity is the driving force for the convective 
motion. A linear relation is more plausible, but at low 
gravity levels we expect the highest deviations, because 
(1) cohesion forces between the particles become impor- 
tant, (2) friction becomes less important due to smaller 
normal forces, and (3) the role of collisions changes due 
to the velocity dependence of the coefficient of restitu- 
tion. Without a numerical model it is, however, difficult 
to predict whether we expect the convection to be en- 
hanced or damped at low gravity levels. 

A linear relation between the convective flux (thus, the 
rise velocity) and gravity was predicted in the model of 
Rajchenbach [12 , which is based on the dilatancy of the 
granular medium and the gradients in density and mo- 
bility. Other models based on thermal convection [14] or 
wall friction [TF describe a 'slight' or square root depen- 
dence on gravity, respectively. With our presented grav- 
ity dependence, we provide a new parameter to verify (or 
falsify) granular-convection models and we can already 



rule out any model with a vanishing or weak dependence 
between convective velocity and gravitational level. We 
should, however, be careful with a too strong statement, 
because we can provide only three data points and are 
therefore unable to discriminate two superimposing ef- 
fects (e.g., convection plus percolation) with a different 
gravity dependence. In fact, the rise velocity for the Lu- 
nar and Martian experiments are in a perfect square-root 
dependence on velocity (dashed line in Fig. HI), and an 
additional superimposed effect could explain the faster 
rise velocity for the terrestrial experiments. 

The strong deviation between our linear fit and the 
rise velocity at the Lunar gravity level inhibits an ex- 
trapolation to much lower gravity levels. A promising 
approach to improve our knowledge on the gravity-level 
dependence of the rise velocity would be experiments un- 
der the enhanced acceleration conditions of a centrifuge, 
which could establish a reliable scaling of the convec- 
tion velocity over orders of magnitude. In the meantime, 
physical and numerical models for granular convection 
should be reconsidered to explain a (roughly) linear de- 
pendence on the ambient gravity. 

We thank Pilipp Reiss, Ralf Pursche, and Martin Rott 
for sharing their parabolic flight with us and for help- 
ing us with the preparation and performance of the ex- 
periments. The parabolic flight was kindly provided by 
ESA and the experimental hardware was funded by DLR 
through grant 50WM1236. C.G. is grateful to the Japan 
Society for the Promotion of Science for his funding. 



[1] A. Rosato, K. J. Strandburg, F. Prinz, and R. H. Swend- 

sen, Phys. Rev. Lett. 58, 1038 (1987). 
[2] J. B. Knight, H. M. Jaeger, and S. R. Nagel, Phys. Rev. 

Lett. 70, 3728 (1993). 
[3] J. B. Knight, E. E. Ehrichs, V. Y. Kuperman, J. K. Flint, 

H. M. Jaeger, and S. R. Nagel, Phys. Rev. E 54, 5726 

(1996). 
[4] J. Duran, T. Mazozi, E. Clement, and J. Rajchenbach, 

Phys. Rev. E 50, 5138 (1994). 
[5] L. Vanel, A. D. Rosato, and R. N. Dave, Phys. Rev. Lett. 

78, 1255 (1997). 
[6] A. Garcimarti'n, D. Maza, J. L. Ilquimiche, and 

L Zuriguel, Phys. Rev. E 65, 031303 (2002). 
[7] W. Cooke, S. Warr, J. M. Huntley, and R. C. Ball, Phys. 

Rev. E 53, 2812 (1996). 
[8] M. E. Mobius, X. Cheng, P. Eshuis, G. S. Karczmar, 

S. R. Nagel, and H. M. Jaeger, Phys. Rev. E 72, 011304 

(2005). 
[9] P. Hejmady, R. Bandyopadhyay, S. Sabhapandit, and 

A. Dhar, Phys. Rev. E 86, 050301 (2012). 
[10] K. Ahmad and I. Smalley, Powder Technology 8, 69 

(1973), ISSN 0032-5910. 
[11] G. M. Rodriguez-Liiian and Y. Nahmad-Molinari, Phys. 

Rev. E 73, 011302 (2006). 
[12] J. Rajchenbach, EPL (Europhysics Letters) 16, 149 

(1991). 



[13] E. L. Grossman, Phys. Rev. E 56, 3290 (1997). 

[14] R. Ramirez, D. Risso, and P. Cordero, Phys. Rev. Lett. 
85, 1230 (2000). 

[15] S. Klongboonjit and C. S. Campbell, Physics of Fluids 
20, 103301 (2008). 

[16] E. Asphaug, P. J. King, M. R. Swift, and M. R. Merri- 
field, in Lunar and Planetary Institute Science Confer- 
ence Abstracts (2001), vol. 32 of Lunar and Planetary 
Institute Science Conference Abstracts, p. 1708. 

[17] H. Miyamoto, H. Yano, D. J. Scheeres, S. Abe, 
O. Barnouin-Jha, A. F. Cheng, H. Demura, R. W. 
Gaskell, N. Hirata, M. Ishiguro, et al.. Science 316, 1011 
(2007). 

[18] B. Gundlach and J. Blum, Icarus 223, 479 (2013). 

[19] The rise time, which determines the rise velocity, was 
calculated by adding up the time intervals At of the air- 
craft's acceleration sensor during the time of excitation 
and individually weight them with [anom/t[,cai(i)]~^'^. 
Here, Onom = flcxc/Samb — 1 = F — 1 is the desired ex- 
citation acceleration and aroai(i) = acxc/<?(t) — 1 is the 
actual excitation relative to the aircraft gravity g(t) at 
time t. The slope of -1.3 is the inverse slope of Eq.flland 
was confirmed in an independent analysis of the rise time 
and excitation acceleration. This correction changes the 
rise times and thus velocities by approximately 10%. 



